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■ Abstract 

wo: 

In this paper, we present extraordinary algebraic and geometrical structures for the 
| Hunter-Saxton equation: infinitely many commuting and non-commuting x, i-independent 

higher order symmetries and conserved densities. Using a recursive relation, we explicitly 
generate infinitely many higher order conserved densities dependent on arbitrary parameters. 
We find three Nijenhuis recursion operators resulting from Hamiltonian pairs, of which two 
are new. They generate three hierarchies of commuting local symmetries. Finally, we give a 
£/3 , local recursion operator depending on an arbitrary parameter. 

As a by-product, we classify all anti-symmetric operators of a definite form that are 
compatible with the Hamiltonian operator D~ x . 
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1 Introduction 

> 
in 

OO . The Hunter-Saxton (HS) equation 

^ ! u x t = 2uu xx + u 2 (1) 

was proposed by Hunter and Saxton as an asymptotic model for the propagation of weakly 
nonlinear unidirectional waves [7]. Its integrability was proved by Hunter and Zheng [8] by 
studying the nonlocal evolution equation 



u t = 2uu x - D x 1 u x , (2) 
where D~ x is the inverse of the total derivative D x . Indeed, equation (2) can be written 

ut = D~ x 5 u {-uu 2 x ) = (u x D~ 2 - D~ 2 u x ) 5 u (-u 2 x ), (3) 

where 6 U is variational derivative with respect to the dependent variable u. Operators D~ x 
and u x D~ 2 — D~ 2 u x are Hamiltonian and form a Hamiltonian pair. This leads to a recursion 
operator of the Hunter-Saxton equation, 

= {u x D- 2 - D~ 2 u x )D x . (4) 

The x-derivative of Hunter-Saxton equation (1), that is 

u xx t — 2,wu xxx + Au x u xx (5) 

is closely related to the Camassa-Holm equation 
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We often write the Camassa-Holm equation in the form: 

mt = um x + 2mu x , m = u xx — u . 
Equation (5) corresponds to TT\i — "^xx 

under the time scaling transformation t >-> 2t. As pointed 
in [8], the bi-Hamiltonian structure (3) and the Lax pair of the Hunter-Saxton equation can be 
obtained from the corresponding known structures of the Camassa-Holm equation. Geometri- 
cally, the Hunter-Saxton equation (1) describes geodesic flow associated to the right-invariant 
metrics on a homogeneous space [9]. It is a particular case of the Euler-Poincare equation on 
the diffeomorphisms in one spatial dimension [6]. 

Recently, equation (2) was proposed as a model to describe shortwave perturbation in a relaxing 
one-dimensional medium. Its integrability was studied by introducing v = D" 1 ^. This leads 
to the hydrodynamic system 

ut = 2uu x — v, vt = 2uv x , 

which was called an integrable regularization of equation (2). We refer to [5, 13] and the 
references in them for the more details in this aspect. 

In this paper, we look at the Hunter-Saxton equation (1) in its own right instead of the traditional 
approach of relating it to the Korteweg-de Vries equation and the Camassa-Holm equation [9]. 
Quite surprisingly, we find it possesses extraordinarily rich algebraic structures: possessing 
infinitely many commutative and non-commutative x, t-independent higher order symmetries 
and conservation laws. Besides the bi-Hamiltonian structure in (3), we found another two bi- 
Hamiltonian structures. Using these operators, we can obtain local recursion operators for the 
Hunter-Saxton equation (1). 

We note that equation (5) can be linearised by the transformation [1] 

v = — , y = u x . (6) 

The linearized equation is 

vt = -y 2 v y - 3yv . (7) 

However, we didn't find a direct way to produce the results in this paper via the linearisation. 
We also note that the calculation in this paper is purely algebraic. We do not justify it in 
analytical sense. 

The arrangement of the paper is as follows: In Section 2, we define the required concepts such 
as Hamiltonian, symplectic and Nijenhuis operators, symmetries, cosymmetries, conservation 
laws and recursion operators for evolution equations in the context of the variational complex. 
We then devote the rest of the paper to the study of symmetries and conservation laws of the 
Hunter-Saxton equation (2). In Section 3, we present a recursive relation to generate infinitely 
many conserved densities (see Theorem 2). We compute their Poisson brackets with respect 
to the Hamiltonian operator D x l and obtain only three commuting pairs among conserved 
densities and T^'^ defined in Theorem 2. In Section 4, we find three recursion operators 
corresponding to three commuting pairs obtained in the previous section (see Theorem 3). To 
prove that these operators are Nijenhuis, we classify all anti-symmetric operators of the form 

U = f(u x ,u xx )D x + D x f + g^.u^u^D^h^^u^u^) + hD x l g, / ^ 0, 

which are compatible to Hamiltonian operator D~ l . Here / is a smooth function of u x and u xx ; 
g and h are smooth functions of u x , u xx and u xxx . We list all five cases that arise (see Theorem 
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4). From one of the Nijenhuis recursion operators, we construct a parameter-dependent local 
recursion operator, which is no longer Nijenhuis (see Corollary 2). Finally, we complete the 
paper with some discussion in Section 5. 

2 Definitions 

In this section, we sketch the basic definitions of Hamiltonian, symplectic and Nijenhuis op- 
erators following [4, 2, 17]. In the context of the variational complex we also define the some 
concepts for evolution equations such as symmetries, cosymmetries, conservation laws and re- 
cursion operators, which also serves to fix our notation. 

2.1 Complex of variational calculus 

Let x,t be the independent variables and u be a (vector- valued) dependent variable. All smooth 
functions depending on u and x-derivatives of u up to some finite, but unspecified order form a 
differential ring A with total x-derivation 

°° d 

D x = ^2 Uk + l T> ' Where U k = d x U - 

k=0 

The highest order of x-derivative we call the order of a given function. For any element g G A, 
we define an equivalence class (or a functional) Jg by saying that g and h are equivalent denoted 
g = h if and only if g — h G ImD x . Without causing confusion we sometimes write g instead of 
J g . The space of functionals, denoted by A', does not inherit the ring structure from A. 

The derivations on the ring ^A. commuting with D x are known as evolutionary vector fields. 
They are of the form dp = Y^T=oD k P-^- Let f) denote the space of all such P. The natural 
commutator of derivations leads to the Lie bracket on f) , that is 

[P, Q] = D Q [P] - D P [Q], P,Qe V , (8) 

where Dq = J2ilo f^^z ^ s tne Frechet derivative of Q. 

The action of any element P € I) on J g <G A' can be defined as 

Pjg = jd P (g) = j IXo = J D g [P\. (9) 

This action is a representation of the Lie algebra f) . We build up a Lie algebra complex associated 
to it. This complex is called the complex of variational calculus. Here we give the first few steps. 

We denote the space of functional n- forms by Q n starting with Sl° = A'. We now consider 
the space fi 1 . For any vertical 1-form on the ring A, i.e., u = YlkLo^duk, there is a natural 
non-degenerate pairing with an element Pel): 

< w, P >= I IXo h k D k x P = f (EZo(-Dx) k h k ) P . (10) 

Thus any element of Q 1 is completely defined by £ = J2T=o(~ D x ) k h k . 

The pairing (10) allows us to give the definition of (formal) adjoint operators to linear (pseudo)- 
differential operators [12]. 
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Definition 1. Given a linear operator S : f) — > fi 1 , we call the operator S* : f) — > S7 1 i/ie adjoint 
operator of S if < SPi, P 2 >=< S*P 2 , P\ >, where Pi e t) for i = 1, 2. 

Similarly, we can define the adjoint operator for an operator mapping from 1 to f), from f) to 
f) or from S7 1 to fi 1 . 

The variational derivative of each functional g G A' denoted by 5 u g G il 1 is defined so that 

< P >= (dfg)(P) =< E2Lo(--D*)*^. ^ > - (11) 

where d : il n — > f2 n+1 is a coboundary operator. Due to the non-degeneracy of the pairing (10), 
we have 

^ = £(-^)'^<eo'. 

fc=0 

In the literature one often uses E referring to the Euler operator instead of 8 U . 

For any £ G fi 1 , by direct calculation we obtain d£ = — D£. We say that the 1-form £ is 
closed if d£ = 0. 

Finally, we give the formulas of Lie derivatives along any K G i) using Frechet derivatives, cf. 
[2] for the details. 

Definition 2. Let Lk denote the Lie derivative along K G f). We have 

L K g = jD g [K] for g G A'; 
L K h = [K, h] for h G f); 

L K §? = L> 5R [K] -D^K + KZ)^ /or »:fj->fj; 
L K n = D n [K] - D K U - HD* K for U : Q l -> fj; 
L X I = Dxl^] + D^I + XL> X for X : t) -»■ fi 1 . 

In this complex we can identify most of the important concepts in the study of integrablc 
systems such as symmetries, cosymmetries, conservation laws and recursion operators. They 
are all characterised by the vanishing of the Lie derivatives with respect to a given evolution 
equation. This will be discussed further in section 2.3. 

2.2 Symplectic, Hamiltonian and Nijenhuis operators 

Definition 3. A linear operator S : f) — » SI 1 (or Q 1 — >■ t)) is anti- symmetric if S = —S*. 

Given an anti-symmetric operator X : [) — >■ Q 1 , there is an anti-symmetric 2-form associated 
with it. Namely, 

u(P,Q) =<1(P),Q >= - <1(Q),P >= -u(Q,P), P,QGF). (12) 
Here the functional 2-form co has the canonical form [12] 

= ^ JduAldu. (13) 



Definition 4. An operator 1 : f} — > SI 1 is called symplectic if and only if the anti- symmetric 
2-form (13) is closed, i.e., dco = 0. 
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It is useful to know that for any £ 6 Q , if £ is not closed, the operator d£ = Z?£ — Z)| is a 
symplectic operator. This can be used to generate new symplectic operators or to determine 
whether a given operator is symplectic or not. 

For an anti-symmetric operator % : Q, 1 — > f), we define a bracket of two functionals / and g as 

{f, g} H =<Suf,n5 u g> . (14) 

Definition 5. The operator % is Hamiltonian if the bracket defined by (14) is Poisson, that is, 
anti- symmetric and satisfies the Jacobi identity 

{{f,g} n > h }n + {{9,h} n , f} H + {{h,f} n , g} n = 0. 

For the Jacobi identity, there are several equivalent formulas given in [2] (see Theorem 5.1). In 
[12] (see p. 443), it is formulated as the vanishing of the functional tri- vector: 

J 9AD H [H9]A9 = 0. 

We are going to use it to classify all the Hamiltonian operators of a given family of operators in 
Section 4. 

Let Jibe & Hamiltonian operator. The Hamiltonian vector fields and their Hamiltonians possess 
the property [2, 12]: 

n5 u {f,g} H = [H5 u f, H5 u g] . (15) 

The Jacobi identity is a quadratic relation for the operator %. In general, the linear combi- 
nation of two Hamiltonian operators is no longer Hamiltonian. If it is, we say that these two 
Hamiltonian operators form a Hamiltonian pair. Hamiltonian pairs play an important role in 
the theory of integrability. They naturally generate Nijenhuis operators. 

Definition 6. A linear operator -ft : fj — > I) is called a Nijenhuis operator if it satisfies 

[fftP, fftQ] - 3fc[3fcP, Q] - K[P, + 5? 2 [P, Q] = 0, P,Q et). (16) 

For a Hamiltonian pair Hi and H2, if Ti\ is invertible, then operator %2Ti-i l is a Nijenhuis 
operator. 

Using the definition of Lie bracket (8), formula (16) is equivalent to 

Lap» = 3fcL P ». (17) 

The properties of Nijenhuis operators [2] provide us with the explanation how the infinitely 
many commuting symmetries and conservation laws of integrable equations arise. In application, 
there are nonlocal terms in Nijenhuis operators. A lot of work has been done to find sufficient 
conditions for Nijenhuis operators to produce local objects [15, 14, 18]. 

2.3 Symmetries and conserved densities of evolution equations 

To each element K G I), we associate an evolution equation of the form 

u t = K. (18) 
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Strictly speaking, one associates to the evolution equation the derivation 



at ^ o d Uk 

As long as the objects concerned are explicitly time-independent as in this paper, there is no 
difference. We refer to [17, 16] for the case when objects explicitly depend on time t. 

Definition 7. Given an evolution equation (18), when the Lie derivatives of the following vanish 
along K G f) we call: g £ A' a conserved density; h £ f) a symmetry; £ € fi 1 a cosymmetry; 
3? : f) — >• f) a recursion operator; a Hamiltonian operator % : O, 1 — > f) a Hamiltonian operator for 
the equation; a symplectic operator X : f) — > SI 1 a symplectic operator for the equation. 

From the above definitions, we can show that if / <G -4,' is a conserved density of the equation, then 
5 u f is its cosymmetry. Moreover, if % is a Hamiltonian operator and X is a symplectic operator 
of a given equation, then HZ is a recursion operator. The operator H maps cosymmetries to 
symmetries while X maps symmetries to cosymmetries. 

We say that the evolution equation (18) is a Hamiltonian system if for a (pseudo-differential) 
Hamiltonian operator H, there exists a functional / € A', called the Hamiltonian, such that 
% S u f is a symmetry of the equation. Additionally, if for a (pseudo-differential) symplectic 
operator X, which is compatible with Ti, there exists a functional g € A' such that 

Zu t = ZK = 5 u g , 

we say that the evolutionary equation is a (generalised) bi-Hamiltonian system. 



3 Conserved densities of the Hunter-Saxton equation 

In this section, we give the recursive relation to generate infinitely many conserved densities for 
the Hunter-Saxton equation T^ ai ' ' ak \ where A; is a nonnegative integer and ati are parameters. 
In general, these conserved densities are not in involution with respect to the Hamiltonian 
operator D^ 1 . We show that there are only three commuting pairs among and T!f n \ 

Proposition 1. Equation (1) possesses a conserved density of the form T = u\ (jjij satisfying 
D t T = D x (2uT), where a is a constant. 

Proof. According to the definition 7, if T is a conserved density, then its Lie derivative along 
equation (1) vanishes. This is equivalent to D t T G lmD x . We have 

= *,,(*». + „?> + a «i (g- 1 - ?f) = d x ( 2 „» ; (-)") 

and thus we proved the statement. o 

Since the Hunter-Saxton equation is free of any parameter, we get new conservation laws by 
differentiating T with respect to parameter a. Thus we have 
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Corollary 1. Expression S = u\ (j^j ( m ^f) , n GN, is a conserved density of equation (1) 
satisfying D t S = D x (2uS). 

In what follows we show how to build up more conserved densities using the above Proposition 
and Corollary. First we prove the following general result. 

Theorem 1. Assume that F is a conserved density of equation (1) satisfying F t = D x (2uF). If 
for all solutions of equation (1) there exists a function r such that rt = 2ur x , then both Fr a and 

G = F (^r) are conserved densities satisfying D t (Fr a ) = D x (2uFr a ) and D t G = D x (2uG) 
for any constant a. 

Proof. We first check that D t (Fr a ) is in the image of D x . Indeed 

D t (Fr a ) = F t r a + uFr^n = D x (2uF)r a + 2uaFr a ' 1 r x = D x (2uFr a ) . 
Next we show that 

D x r\ a \ ( D x r\ a 1 / D x rt u x tD x r 



n I ! \ i • I -4r I I — ^ -2 



D x r\ a 1 ( D x (2ur x ) (2uu 2 + u\)D x r 

n[ -uT) V^^" 2 5 



n [D x r\ a ( n u 2 D x r\ ( D x r 

= 2au {-uf) y-uj- ^r) = x \u{ 

Using the first part of the proof, we obtain that G is a conserved density for any constant a. o 

We now search for function r such that r t = 2ur x for equation (1). Notice that r satisfies 
r t = 2ur x if and only if lnr satisfies the same relation, that is, (lnr)t = 2u(lnr) x . 

Proposition 2. Assume that both T and Tr are conserved densities of equation (1) satisfying 
T t = D x (2uT) and D t (Tr) = D x (2uTr). Then r t = 2ur x . 

Proof. From the assumption we have 

= D t {Tr)-D x (2uTr)=T t r + Tr t -D x {2uT)r-2uTr x =T(r t -2ur x ). 
This leads to the conclusion that r t = 2ur x . o 
From Proposition 1 and Corollary 1, it follows that r = In ^ satisfies r t = 2ur x and so does 
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Using Theorem 1, we can now obtain the following result: 



Theorem 2. The functionals y( ai ''"' afc ) generated by the recursive relation: 



T = ul 


^0 = 


1 


T[ ai) = T r°\ 


n = 


u{ 


T ,(ai,a 2 ) _ t (qi) a 2 
1 2 — 1 l '2 ' 


T2 = 


-o-Dxn; 

u{ 


rp(ai,— ,ak) rp{ai,— ,afc_i) a k 
1 k — 1 k-l ' k ' 


Tk = 


■\D x r k _ 
u{ 



1; 



(19) 
(20) 
(21) 

(22) 
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are (fc + l) th order conserved densities of equation (1). Moreover, = ; if 1 ' 1 ^ = and 

when ctk ^ — 1 and k > 2, we have 

fc-i 

rp(ai,— ,Q fe ,l) _ y ^-,(ai,--- ,ai-i,ai-l,ai + i+l,— ,afc-i,gfc+l+^,fc-i) 

i=i K 

Proof. We only need to prove the second part of the statement. It is easy to see that 

Tf } = D x r = 0, T 2 (ai ' 1} = r^D xn = 0. 
We now show formula (23) by direct calculation: 

r £i,..,a fc> l) = Tj («i.-,a fc ) rfc+1 = r m r a 2 . . . 
fc-1 

\ - Qj ai g i _ 1 a fc _i afc +l 

— Zj a I J 1 r i— 1 r i IMz^Wi r k-l r k 



U l r l r i-l r i „,2 'fc-1 'fc 



Thus we obtain the formula using the relation = and the recursive relation for o 

For each A; G N, the above T^ ai ''"' ak ^ is a (fc + l) th order conserved density with fc parameters. 
In a similar way, we can build up more conserved densities by adding more logarithms in front 
of rj. For instance, both T" 1 (lnri f 1 and If 1 (Innf 1 (In In ri) 71 are conserved densities. 

Here all T/% are local. By no means are the conserved densities we constructed above complete. 
For example, we haven't include the conserved density that is, 

1 2?/ 
D t (-) = D x (—-3x). 

U\ Ul 

Besides, there are also nonlocal conserved densities. In [8], the authors listed some of them 
generated by the recursion operator (4). For example, the conserved density 2u 2 u\ + (Z)" 1 ^) 2 
depends on the same nonlocal term D" 1 ^ 2 as in equation (2). 

We can define the Poisson bracket of any two conserved densities with respect to a Hamiltonian 
operator D^ 1 according to formula (14): 

{•/•;- - --' . if 1> - A) } x =< 5 u (7f 1, '"' a<) ), D-H u {T^'-' Pi) ) > . (24) 
It is clear that 

{t , T^'-' a ^] D i =<-2u 2 , D- 1 S u (T^'-' a ^)>=<2u 1 , «5 u (rf 1 '-'^ ) ) > 

=< 2u u D*, ait ..., aj) (l)) >= [ 21) ■ = . (25) 

7 ^ 
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(26) 



With the help of the computer algebra system Maple , we obtain that 

| r («) > jW) }^ = l a/3(Q( _ 1)G g _ 1)(a _ /3) r (« + /3-5,3). 
| T H ; = l a7(a _ 1)(7 _ 1)(2 _ 7) T (« + /3- 2 , 7 -3,3) 

+ I a7 ( a _ 1)( 7 _ l)( 5a _ 3/3 _ 10) r^-3,7-1,2) 

+a(a - 1)( 7 - 1)(6 7 + 7a 2 + /3 - a/3 7 + 2/3 7 - /3 2 - 5 7 a) r^-^+i.i) 
+a/3(a - l)(a - 2)(1 - /3)( 7 - 1) T 2 (a+,9 " 5 ' 7+3) . (27) 



This implies that nontrivial (i.e. a(a — 1) 7^ 0) and T^ p> do not commute unless a = f3. 

Now we look at when a nontrivial commutes with a nontrivial T 2 ^' 7 ' ) . From Theorem 2, we 

know that T 3 (ai ' a2,1) = ^^r 2 (Ql_1 ' Q2+2) . Hence, the right-hand side of formula (27) vanishes if 
and only if 

f 2-7 = 0; 

5a - 3/3 - 10 = 0; 

-^^(67 + 7" 2 + /? " a/37 + 2/?7 - /? 2 - 57a) + (a - 2)0(1 - /3) = . 
Solving this algebraic system, we obtain the following three solutions 
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a = 2 
/? = -§ 
I 7 = 2 



a = — 1 
= -5 
I 7 = 2 




These lead to the following statement. 

Proposition 3. Assume that a(a — 1)7(7 ~~ 1) / 0- 

among t[°^ and T 2 ^' 7 ' ) , namely, 



There are only three commuting pairs 



2\2 



U 



1/2 (uiuz - 4n 2 ) 

2 5/2 
«1«2 



) K 

^2 



u 2 («i«3 — 4n 2 ) 2 



6 " 
1 



u 



. 



DZ. 



Note that the Poisson bracket of any pair of conserved densities is again a conserved density 
since 

It is surprising to see the results are linear combinations of the conserved densities listed in 
Theorem 2 since they are by no means complete. It seems that these conserved densities are 
closed under the defined Poisson bracket (24). 

In next section, we show that the above three commuting pairs lie in three different commuting 
hierarchies generated by three Nijenhuis recursion operators. The natural question is whether 
there are more commuting pairs if we compute the Poisson bracket between higher order con- 
served densities. In Appendix A, we include the formula for |t i ^' ) , T^' 1 '^ j _ i . Based on it, 

we find three commuting pairs between conserved densities of second order and those of fourth 
order, generated by the same three Nijenhuis operators. We conjecture that there are only three 
commuting hierarchies starting with T^ a \ By computing the Poisson bracket between j 1 ( ai ' a2 ' ) 
and j 1 ^ 1 '^ 2 '^ 3 ^ we discover another commuting pair. It is listed in Appendix A. However, we 
have not found the corresponding Nijenhuis recursion operator. 
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4 Symmetries and Recursion operators 



We know from (3) that the operator D~ x is a Hamiltonian operator for equation (2), mapping 
cosymmetries (the variational derivatives of conserved densities) to symmetries. Thus, we can 
produce infinitely many symmetries from the conserved densities listed in Theorem 2. Using the 
property of Hamiltonian operators (15), we have 

D-'Su {•/;" -••'". T? 1 '-'™}^ = [LT^Crf 1 ''"' Qi) ), D-'duiT^ 1 '-'^)}. (28) 

From the results in the previous section, we know some symmetries are commuting and some 
are not, cf. formula (26) and (27). In this section, we will present some recursion operators 
to generate infinitely many commuting and noncommuting symmetries of the Hunter-Saxton 
equation (2). 

4.1 Nijenhuis recursion operators and commuting symmetries 

The Nijenhuis recursion operators [11] are used to generate infinitely many local commuting sym- 
metries for integrable equations. In this section, we present three Nijenhuis recursion operators 
for the Hunter-Saxton equation (2), which correspond to three commuting pairs in Proposition 
3. We then prove that they generate infinitely many local symmetries. 

Theorem 3. The following three operators: 

»i = (2u 2 l D x + 2D x u 2 1 - u 3 u 2 3/2 D x 1 u 3 u 2 3/2 ^j D x (29) 

sft 2 = (u\u 2 2 D X + D x u\u 2 2 - 8mD- V) D x (30) 
sft 3 = ( U ^D X +D x u^- 4(u 3 ur 6 - 3ulu^ 1 )D-\i l - 4 Ul D-\u 3 u^ 6 - ?>u\u{ 7 ))D X (31) 

are all recursion operators of equation (2). 

Proof. According to Definition 7, we check whether Li^iti vanishes for each i = 1,2,3, where 
K = 2uu\ — D~ x u\. Since the calculation is similar, we only work it out for i = 2. To simplify 
the computation, we introduce some notation: s = u\u 2 2 and Sj = D x s. So 

= 2{uD x + Ul - D- l uiD x ){sD x + D x s - 8u l D- 1 u l )D x 

-2{sD x + D x s - 8u 1 D x 1 u 1 )D x {uD x + Ul - D x x u x D, x ) 
= 2{uD x + D~ l u 2 )(sD x + D x s - 8uiD x 1 ui)D x 

-2{sD x + D x s - %uiD x l ui){uD x + m + u 2 D x 1 )D x 
= ((2nsi - 4sui)D x + D x (2usi - Asm) - 8(2uu 2 + u\)D~ l ui + D^f^ul - 4su 3 - 2siu 2 ) 

-8u 1 D- 1 (2uu 2 + u\) + (8u\ - Asu 3 - 2s l u 2 )D x : 1 ) D x 
= (A{2uu\u 2 1 — uufu 2 3 u 3 — u\u 2 2 )D x + AD x {2uu\u 2 ~ 1 — uufu 2 3 u 3 — u\u 2 2 ) 

-8(2uu 2 + Ui)D~V - 8u 1 D x 1 (2uu 2 + u\)) D x , 

which equals to D^ 2 [K]. Therefore, we have Lj^fi 2 = and thus the statement is proved. o 

The recursion operator (4) given in [8] is the inverse operator of Indeed, we can prove the 
following statement: 
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Proposition 4. ftj = 85ft- 1 = SD^^D- 2 - D^m)' 1 . 

Proof. Note that u\D~ 2 — D~ 2 u\ = D~ 1 (u2D~ 1 +D~ 1 U2)D~ 1 . To prove the statement, we only 
need to show that 8 (u2D~ l + D X 1 U2) 1 = 2u 2 1 D x + 2D x u 2 x — u 3 u 2 3 ^ 2 D x 1 u 3 u 2 3 ^ 2 ■ Indeed, 

(2u 2 1 D x + 2D xU2 ~ 1 - u 3 u~ 3/2 D'^u'^^D- 1 + Z^V) 
= 8+2u 2 1 u 3 D x 1 -2u 2 2 u 3 D x l U2-2u 3 u 2 3/2 D x l {u 2 ,2 ) x D x 1 + 2n 3 n^ 3/2 Z) a : 1 (n^ 1/2 ) a; D- 1 n2 

= 8 + 2u 2 1 u 3 D- 1 - 2u 2 2 u 3 D- 1 u 2 - 2u 3 u 2 1 D~ 1 + 2u 3 u 2 3/2 D^uJ 2 

o -2 r>-l o - 3 / 2 n-1 V 2 

+2u 3 u 2 D x U2 — 2u 3 u 2 D x u 2 

= 8, 

and this leads to the statement. o 

In the paper [8] , the authors proved that the recursion operator (4) is the ratio of a Hamiltonian 
pair. Hence it is a Nijenhuis operator and so is its inverse operator [2]. We now prove that 
the other two recursion operators 3?2 and $t 3 are also the ratio of Hamiltonian pairs. To do so, 
we classify all anti-symmetric operators of the form 

H = f(u 1 ,u 2 )D x + D x f + g(u 1 ,u 2 ,u 3 )D x 1 h(u 1 ,u 2 ,u 3 ) + hD^g, / / 0, (32) 

which are compatible with the Hamiltonian operator D~ x . Here / is a smooth function of u\ 
and U2] g and h are smooth functions of u\, U2 and u 3 . We give the result below and the proof 
in Appendix B. 

Theorem 4. If an anti- symmetric operator (32) forms a Hamiltonian pair with the Hamiltonian 
operator D^ 1 , then the smooth functions f , g and h are one of the following five cases. Here 
a, b, c, d are smooth functions of u\ and qCC, i = 1, 2, • • • , 5 are constants. 



I. 



II. 



III. 



IV. 



— . «.(•«,) / 



f a(ui) 2 u 2 

h = b(u\) 
9 = c(ui) 

f = 1 

•> (<j(ui)«2+Ci) 2 

h = c 2 u x + c 3 
g = C4U1 + C5 

f = 4 

(a(ui)u2+b(ui))' 1 ' 

I = C\U\ + c 2 

16tt3 



(33) 



(34) 



o(«i) + 



- ■ 1 g_9_ ( b+2au 2 

y (a(ui)u2+b(ui))' A dui \a(u\) 2 (a(ui)u2+b(u\)) 2 

b(ui) = c\u\ + 2c 2 ui + c 3 

b'c 4 +(nife / -fe)c 5 16a'fc+8afe' 



)+c(«i) (35) 



c(iti) 

/ = 
ft = 



b(«i) 2 



o(ui) 3 6(ui) 2 



5 = " 
6(«i) 
c(ui) 



C\U\ + c 2 
16m, 



. 2ib'u 2 

H 



b( Ul ?^^s^ + « 
= c\u\ + 2c 2 «i + c 3 

fe'c4+(uife'— b)c5 



(36) 
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a(ui)u|+6(ui)u2+c(«i) 



~~ (a(ui)u|+6(ui)M 2 +c(«i)) 3 / 2 2cic 2 &uT \ (a(ui)u|+6(ui)M 2 +c(«i)) 1 / 2 

5 = c 2 /i - 2c 2 cZ(ui) 

= c 3 ni + c 4 
4ac - 6 2 = 4cfc 2 




V. 



< 



(37) 



Example 1. Two special cases from case II (cf. (34)) lead to Hamiltonian pairs 



D x + 2u 1 D x x + 2D x l Ul + XD 



and 



D x + 2u x D- x u x + AD' 1 . 



These are the hi- Hamiltonian structures for potential Korteweg-de Vries equation and potential 
modified Korteweg-de Vries equation respectively since 



The recursion operators (30) and (31) in Theorem 3 are the ratios of Hamiltonian pairs listed 
in Theorem 4. It is obvious that the recursion operator (30) can be derived from Case I (33). 

For Case IV, if we take b = 2u\, c{u\) = and h = 2u±, then we get g = —2ui 6 us + 6ui 7 ul. 
This leads to a Hamiltonian pair 

u^D x +D x u^ 4 - 4(u 3 Ui 6 - 3ulu^ 7 )D x 1 ui- AuiD x Xuzu^ - 3u 2 2 u~{ 7 ) + AD" 1 , 

from which we obtain the recursion operator 5ft 3 (31) in Theorem 3. We can also directly get the 
recursion operator (29) from Case V by taking a{u\) = c(u\) = d(u\) = and b{u\) = ^. Since 
all three recursion operators can be obtained from Hamiltonian pairs, they are all Nijenhuis 
operators. 

Theorem 5. The recursion operators (29)-(31) are Nijenhuis operators. 

We can also view the recursion operators (29)-(31) as the products of Hamiltonian and sym- 
plectic operators. In a recent paper [18], we proved that for Nijenhuis operators that are the 
products of weakly nonlocal Hamiltonian and symplectic operators [10], hierarchies of commut- 
ing local symmetries and conserved densities in involution can be generated under some easily 
verified conditions. To be self-contained, we restate the result in [18] valid for the operators in 
this paper: 

Consider a Hamiltonian operator H of the form (32) and a symplectic operator 
T = D x such that "HZ is a Nijenhuis operator. Assume that L g T = L g Ti = L^T = 
LhH = 0. If there exists a closed 1-form £ satisfying L g ^ = = = such 

that XHi is closed, then all {X%y^ are closed 1-forms and commute for 



and 




i = 0,l,2,---. 
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We denote the recursion operators (29)-(31) in Theorem 3 as 9?j = HiD x , i = 1,2, 3. Checking 
the conditions in this statement for all ift^, we can prove the following theorem. 

Theorem 6. Let £i = 5 u (u^ 2 ), £2 = ^u{^) and £3 = <5 u (^lr). Then all are closed 

1-forms and for each fixed i € {1, 2, 3}, all H^Jt* 1 ^ commute for j = 0, 1, 2, • • • . 

Proof. We write out the proof for the recursion operator 3?2- The proof for operators 3?i and 
is similar, and will not repeat it. Since £2 is the variational derivative of it is clearly a 
closed 1-form. In this case, we have g = h = u\ and trivially L Ul D x = L Ul %2 = L Ul ^2 = 0. We 
only need to show that L^ 2 ^2 = and D x %2^,2 is closed. Note that 

i2 = Hfy _ 3^1 + 12^ - 30.?; 

u 2 U 2 U 2 

U2i2 = 2^ - 20^i + 40^i + 30^1 - 120^1 + 220^ - 320^; 

«2 11% u\ U' 2 u\ U2 U2 



and D x %2^,2 = —5 u { U ^ UlUi h ), which implies that D x %2£,2 is closed. Using the computer 
algebra system Maple , we can check L^ 2 ^2 = D^ 2 [H^2] + ^^^(^2) = 0. Thus we prove the 
statement for 3?2- o 

We have proved that the recursion operators (29)-(31) are Nijenhuis. Using them, we generate 
three hierarchies of commuting symmetries. However, the elements in the different hierarchies 
do not commute. 

Note that the Lie derivative is a derivation. Thus the products and additions of recursion 
operators are recursion operators. For instance, operators K1K2 and [Sfi,^] are also recursion 
operators of the Hunter-Saxton equation (2). In general, they are no longer Nijenhuis (cf. (16)) 
and do not generate hierarchies of local symmetries. 



4.2 Recursion operators and noncommuting symmetries 

In this section, we first prove that the adjoint operator of the recursion operator K2 (30) acting on 
any cosymmetry 5 u (Tj ai ' ' a ^) produces local cosymmetries. Since the resulting cosymmetries 
are not closed, we can construct symplectic operators, which depend on parameters. Further, 
this leads to parameter-dependent recursion operators. 

Proposition 5. Starting from any symm etry Q = D x 1 5 u (Tj ai '"' ' a ^) ; where T^ ai ''"' aj) ) are 
listed in Theorem 2, all Q k = St^Qo are local symmetries k = 0, 1,2, Equivalently, all 
$t,2 k $u(Tj ai ' ' aj ' 1 ) are local cosymmetries, where ffi, ^ s ^ e adjoint of operator 0^2 • 

Proof. Note that 5?2 = {u\u2~ 2 D x + D x u\u2~ 2 — SuiD^ui) D x . To prove that Q k are local, we 
only need to show that U2Qk is in the image of D x . We prove the statement by induction. 

We know that from (25) the Poisson bracket of To and Tj" 1 ' vanishes. This implies that 
U2Q0 G lmD x . Assume that U2Qk-i £ lraD x . We now show that U2Qk € lmD x . Indeed, 

U2Qk = u 2 ^2Qk-i = U2 (ufu2~ 2 D x + D x u\ii2 2 - SuiD' 1 ^) D x Q k _i 

= 2ufu2 1 DlQ k -i - 2ufu 2 2 u 3 D x Q k ^ 1 + <iu\D x Q k _i - Su\u 2 Q k ~i + 8uiu 2 D~ 1 u 2 Qk-i 
= D x (2uju2 1 D x Q k _ 1 - Au\Q k -i + Au\D x l U2Q k -\) € lmD x . 
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Thus all Qk are local. Note that 

u 1 S% k 6 v (Tj ai ''"' a ' ) ) = Ul D x n 2 ■ ■ ■ D x H 2 6 u (T^ ai ' - ^ ] ) = u x D x Q k = D x ( Ul Q k ) - u 2 Q k • 

We have just proved that u 2 Q k € lmD x . Thus we have uiS u (Tj ai ' ' aj ^) € lmD x . Hence all 
$t6 u (T$ aw "' aj) ) are local. o 

In general, the symmetry Qq defined in Proposition 5 is not a symmetry of operator $t 2 , that 
is, Lq 5?2 7^ 0. Therefore, [QojQi] 7^ 0. Although $l 2 is a Nijenhuis operator, the generated 
symmetries do not commute. Furthermore, $Q8 u (Tj ai ' ' a ^) are no longer closed. 

Let us look at a simple case when j = 1. The corresponding cosymmetry is 

5 u T[ a) = 5 u {u{- Aa u%) = (q - l)D x (auj- 4a u%- 2 u 3 + (2 - 4a)u\- 4a u%) . 

Notice that 

5 u {u^ 2 u 2 ln(u^ A u 2 ) = D x (u^u^us - 2u± 3 u 2 ) . 

Without losing generality, we let tft^ act on the cosymmetry 

= D x (auj-^u^us + (2 - Aa)u\-^u% ) , 

where a <G C and £W is closed. We have 

sft^ 1 ) = 2aD x (ul- 4a u%- 4 u 5 + (3a - 7)u\- ia u^u z u A + (8 - 12q)^- 4q uJTV 
+(a - 4) (a - 2)u?- 4a u£- 6 ujj - 4(3a - 2)(a - 2)uf- 4a uJ- 4 ui 
+2(4a - l)(6a - 5)u 4 - 4Q n°- 2 n 3 - 16(4a - l)(a - l)n?" 4a uf) , 

which is a cosymmetry of equation (2), that is, L^^Si^ 1 ^) = 0. As we mentioned in Section 2, 
d(3?2^ 1 ^) 7^ is a symplectic operator. We have the following. 

Proposition 6. Operator 1 = D X (D X + ^) ( ^p^ D x + D x tllTi ) (A* - ^)A* is a 
symplectic operator for the Hunter-Saxton equation (2) for all a G C. 

Proof. By direct calculation, we have 

d(^«) = D^ ( i)-D^ (1) 

= (a + 1)D X (D X + — ) 4a „ 5 5 _ a Ar + As 4a-5 5-a Ac ) D * ' ( 38 ) 

which is a symplectic operator when a/-l. Since the Lie derivative commutes with d, that is 

L U( d(^ (1) ) = dL Ut (^ (1) ) = 

implying that it is a symplectic operator for the Hunter-Saxton equation (2). When a = — 1, 
we can write 

l = d(D x S). 

Here S = ^ - + f ^ + 16^ - 43^# + 66^ - 96^ is a symmetry of the 

^2 ^2 ^2 ^2 ^2 

Hunter-Saxton equation (2), which can be verified according to Definition 7. By the same reason 
as above, we prove that I is a symplectic operator for the Hunter-Saxton equation. o 
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For formula (38) when a = — 1, we have d^^ 1 ^) = 0, which is consistent with the result in 
Theorem 6. By the Leibniz rule for the Lie derivative, any product of Hamiltonian and symplectic 
operators of the Hunter-Saxton equation is a recursion operator. Using its Hamiltonian operator 
D~ x and a symplectic operator X as in Proposition 6, we obtain the following result: 

Corollary 2. Differential operator (D x + f ) ( fa-. ^ + Ar tllTl ) (Ac - " a 

recursion operator for the Hunter-Saxton equation (2) for a G C. 

Here we can use different Hamiltonian operators in Theorem 3 instead of D~ l . However, the 
other Hamiltonian operators will bring in nonlocal terms in the recursion operators. Further 
study is required to determine whether the resulting operators produce local symmetries or not. 

5 Discussion 

It is well known that the integrable equations possess an infinite number of commuting conserved 
densities and generalised symmetries. In this paper, we present a new feature for the integrable 
Hunter-Saxton equation: infinitely many noncommuting x, i-independent conserved densities 
and symmetries. We found three Nijenhuis recursion operators and a local parameter-dependent 
recursion operator. We believe that there are more Nijenhuis recursion operators related to the 
conserved densities listed in Theorem 2 since we have found a new commuting pair (39) in 
Appendix A. Note that we can define the Poisson bracket (24) with respect to %2 instead of 
D^ 1 . It will be interesting to extend the study in the paper and to see whether there are new 
commuting pairs. 

The conserved densities in Theorem 2 give rise to infinitely many cosymmetries of the Hunter- 
Saxton equation (2), which are closed 1-forms. In Section 4.2, we showed that the results of 
recursion operator 3^2 (30) acting on such closed cosymmetries are no longer closed. Hence, 
we can generate a lot of local symplectic operators such as in Proposition 6. This will lead 
to local recursion operators as in Corollary 2. The immediate questions are: what are the 
relations among such recursion operators? Can we write down neat formulas for them? Using 
the computer algebra system Maple , it is not hard to compute these operators although the 
expressions are huge. The problem is to present them in a compact way, e.g. as the product of 
1 st order differential operators as in Proposition 6. 

We know the set of symmetries is a Lie algebra under Lie bracket (8). For noncommuting 
symmetries of the Hunter-Saxton equation, we can use them to generate higher order symmetries. 
They are different from master symmetries [3], which generate commuting symmetries. 

The fundamental question is: where do such rich structures for the Hunter-Saxton equation 
come from? As mentioned in the beginning of the paper, the equation is linearizable and has 
a trivial dispersion law. However, the transformation (6) is highly nonlocal. We didn't find a 
direct way to produce the present results from the linearized equation (7). It would be very 
helpful for getting a clear and complete picture if we could find the direct link. 

Appendix A 

In this Appendix, we give the Maple result of the commutator of conserved densities between 
T 1 (o) and T { / n '^. We also list the new commuting pair between T 2 {qi ' Q2) and T^ 1,Mz) . 
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Using the computer algebra system Maple , we obtain the following formula: 
{Tf } , TV™)} = l -a{a - 1)^ - 1)(2 - //) T ^~^-^) 

_l a ( a _ i)^ _ i)( 3 /3 - 7a + 14) T4 (°+0-3.7+i,m-2,2) 

+a(a _ 1) 7(7 _ 1 )(M _ l) (/x + 1) ^«+/»-2,7-2^1,l) 

+a(a — l)(2^ 2 /37 — + fi 2 ^ 2 a — Aa^fi 2 + 8// 2 7 + 3ajfi — 27 2 // 2 — 67// 

+// 2 /3 + 6// + 07 + ///3 + 2 7 2 - 6// 2 + 3// 2 a - Q7 2 - 2 7 - 3//a) t 4 {q+/3 ~ 3 ' 7 ' m,1) 
+a 7 (a - 1)(7 - 1)(7 - 2)0* - 1) T 3 ( ^- 2 ^ +3) 

+07(0 — 1)(— 27 2 /i — cry 2 + cry 2 ^ + 27 2 + 37/3// — 67 + 3«7 — 2/?7 + 67// 
-3a 7 // + 2a// -20 + 2/3-4^ + 4) t^+^.t-i,/^) 

+a(a - 1)(72 + 2 7 /?// + 13/3 + 24 7 - 72// - 60a + 4a 2 7 + /3 2 // + 2/3 2 7/ u 
+a 2 7 2 // — 4a 2 7 / u — a 2 7 2 + a/37 2 // + 4a/3// — afij 2 — 8a/3 — 12a 2 // 
-2a^7/x + 60a// - 9/3// + 3/3 2 + 12a 2 - 67 2 - 24 7 // + 2/3 7 2 + 5a7 2 

-2/3 7 2 // + 20a 7 // - 20a 7 + 67V - 5a 7 V) T^ + ^ 1+W) 
+a(a - 1)(240// - 10//a 3 - 8a/3 - 54a 2 - 144 + 24 7 + 9a 2 7 + 10/3 + /3 3 + /3 2 

+156a - 260a// + 2a 2 /3 + 6a 3 - 2a/3 2 + 90//a 2 - 26a 7 - 7a 3 ) T ( a +^- 5 ^+ 3 ^) 
+a(a - l)(a - 2) (a - 3) (a - 4) (5 - a)// jj»+P-^+5,i^-i) 

Combining with formula (27), we can obtain the following result: 

Proposition 7. Assume that a(a — l)//(// — 1) 7^ 0. There are only three commuting pairs 
among and the third order conserved densities, namely, 

j r (l/2) T (-7/2,0,2) _ 35 T (-11/2,4)1 = f T (-l) T (-7,0,2) _ „ r (-9,4)l 

I 1 ' 3 16 2 J D -i I 1 ' 3 2 J^ 1 

= r T (2 )j T (o,o,2)| _ i=0 _ 

Proof. We search for the linear combinations of T^' 1 '^ + a\T^ 1,11 \ which are in involution of 
under the Poisson bracket (24). From Theorem 2, we know that 

rp{vi,V2,l) _ ~ U 1 ^(1^1-1,^2+2), 

3 ^2 + 1 ' 

^(^1^2,^3,1) _ V l ^(^1-1,^ + 1,^3 + 1) _ ^2 ^(^1,^2-1,1/3+2) 

4 ^3 + 1 3 ^3 + 1 3 

We substitute the above formula into jT^, T^' 1 '^ + aiT^ 1 ' 71 ^} ! and write it as the com- 
bination of independent terms. Note that it vanishes if and only if all coefficients of independent 
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terms are equal to zero. We immediately get \i = 2 and 7 = 0. The other conditions are 



7a - 3/3 - 14 = 0; 
/3 = ft + 2; 
7i = 4; 

(-84 + 70a + 3/3 2 - 4a/3 + 5/3 - 14a 2 ) + 2ai( 7 i - 1)(2 - 71) = 0; 
336 - 14a 3 - 8a/3 + 10/3 + /3 3 + /3 2 - 364a + 2a 2 /3 - 2a/3 2 + 126a 2 

+6ai(5a - 3ft - 10) = 0; 
3ai(a - 2)/3i(l - ft) - ±(a + ft - 4) (2(a - 2)(a - 3)(a - 4)(5 - a) 
+3ai(24 + 4a 2 + ft - 4aft + 8ft - /3 2 - 20a)) = . 

Solving this algebraic system, we obtain the following three solutions 



a 



ft - ~T 
71 = 4 



a = — 1 
ai = —7 

ft = -9 
7i = 4 

U = -7 



a = 2 
ai = 
ft = -2 
71 = 4 
I P = 



These correspond to the commuting pairs listed in the statement. o 

Notice that we get the same values of a as in Proposition 3. In section 4.1, we show that there are 
three Nijenhuis recursion operators corresponding to each value of a. These three commuting 
pairs can be directly found from the corresponding recursion operators. For example, 



-25, 



2) 



7 T. 



(-9 



Proposition 7 implies that there are no other commuting pairs between and conserved den- 
sities generated in Theorem 2 of third order. We conjecture that there are only three commuting 
pairs between and conserved densities generated in Theorem 2 of any higher order. 

We compute the commutator between T^ ai ' a2 ^ and j"^ ai ' a2 ' as \ Beside the three pairs directly 
obtained from Proposition 3 and Proposition 7, we also find the following new commuting pair: 



f T (0,l/3) T (0 -7/3,2) 1 j f U1U3 ~ 4ui 



1/3 



lfu4 — UU1U2U3 + 28n^) S 

u^ 3 (uius — 4ti|) 



2\7/3 



= 0. 



(39) 



We have not found the corresponding Nijenhuis recursion operator for this new case as we did 
in section 4.1. 



Appendix B: Proof of Theorem 4 

In this Appendix, we give the proof of Theorem 4. Let us first introduce some notation: 

fi = Dif, f {j) = ^ U0)-i = D-\fe). 

The same notation is also used for g and h. 

Proof. To prove the statement, we check when the operator 

H = f(ui,u 2 )D x + D x f + g(u 1 ,u 2 ,u 3 )D~ 1 h(ui 1 u 2l u 3 ) + hD^g + \D~ l 
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is Hamiltonian for arbitrary constant A € C. From Chapter 7 in [12], We know % is Hamiltonian 
if and only if 

PrVw(0 w ) = O, where Q H = ± J (9 A H9) dx 

is the associated bi- vector of T~L. 
First we have 

U(6) = 2/0! + hO + gihB)-! + %0)_i + A0_i; 

D X (U{9)) = 2f9 2 + 3/i0i + f 2 + 2gh9 + gi(M)-i + + A0; 

D 2 (H(0)) = 2/0 3 + 5/i0 2 + 4/ 2 0! + / 3 + 2gh6 l + 3( 5 /i)i0 + </ 2 (/i0)_i + h 2 {g9)- l + A0 i; 

D X (H(8)) = 2/04 + 7/i0 3 + 9/ 2 2 + 5/ 3 0i + / 4 + 2gh6 2 + 5(^)i0i + l(gh) 2 9 - 2g 1 h 1 8 

+gs(he)- 1 + h 3 (ge)- 1 + \e 2 . 

Substituting them into Pr Vue{®u)i this leads to 



Pr W(e) e w = J (0 A £>/[W(0)] A 0i + A £> ff [W(0)] A (/i0)_i + A D^(0)] A (<?0)_i) dx 
= y A (f^D x {U{9)) + /( 2 )^(^(0))) A x dx 

+ J 9A(gWD x (H(9)) + gWDl(H(9)) + gWD 3 x (H(e))) A(h8)^dx 
+ J A (h m D x {U{6)) + D 2 X (U(9)) + h^Dl{U{9))) A {g9)- 1 dx . 
It needs to vanish for all A. So the coefficient of A should be zero, that is 

= J A (# {2) 0i + 5 {3) 2 ) A (/i0)_idz + j 6 A (/i (2) 0i + /i {3) 2 ) A (£r0)-idx 
= J (0 A (gW - g ^)e l A (M)-! + A (h^ - / l f))0 1 A (<0)_i) dx . 

This leads to = hf ] and g^ = gf\ The implies that 

h(ui,U2,u 3 ) = Ph(ui,u 2 )u 3 + qh(ui,u 2 ), g(ui,u 2l u 3 ) = p g (ui, u 2 )u 3 + q g (u 1 ,u 2 ) (40) 
and 

?i 2) =^2, ^ 2) =#n 2 . (41) 
The rest of the terms, i.e., the terms without A in PtVho(@h) should also vanish. 

= y'/( 1 )0A(2/02+ 5 i(/i0)_ 1 + / l i( 5 0)_i) A0idx 

+ J f®0 A (2/0 3 + 5/i0 2 + g 2 {h9)^ + h 2 (g9)_ 1 ) A 0idx 

+ J gW0 A (2/02 + 3/i0i + A (/i0)_idx 

+ y 5 ( 2 )0 A (2/0 3 + 5/i0 2 + 4/ 2 0i + 2gh9 1 + h 2 (g9)_ 1 ) A (he^dx 

+ Jg^e A (2/04 + 7/i0 3 + 9/202 + 5/ 3 0i + 2 5 /i0 2 + 5( 5 /i)i0i + /i 3 (^)-i) A (/i0)_ x dx 
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+ J hVO A (2/0 2 + 3/i0i + 9i(h6)-i) A {g0)-idx 
+ J h®9 A (2/0 3 + 5/i0 2 + 4/ 2 0i + 2gh9 l + g 2 {h6)^) A (<70)_idr 
+ Jh®6 A (2/04 + 7/i0 3 + 9/ 2 2 + 5/ 3 0i + 2gh9 2 + 5(^)i0i + g 3 (h6)-i) A (s0)-idz 
= y (2D x (f (2) f) ~ 5/i/ (2) - 2//« + 2//u,< 3 > + 2/^ {3) ) A 0i A 2 dx 

+ J (-i53(2/^ 3 )) + ^(2/ 5 ( 2 ) + 7/i 5 (3)) + D x (-2fgW _ 5 / l5 ( 2 ) - 9/ 2 2 (3) _ 2 ghg^) 

+4f 2 g {2) + 2^^ + 5/ 35 {3) + 5(^)i5 (3) +3/i5 (1) -5i/ {1) -52/ (2) ) A 0i A (h6)^dx 
+ J (-^(2/^(3)) + l£(2/fc< 2 > + 7M(3)) + D X (_2/*W - S/i^ 2 ) - 9/2^3) _ 2 ghhW) 

+4/ 2 /i (2) + 2 5 W 2 ) + 5/ 3 /i (3) + 5( 5 /i)i/i {3) + 3/i/i (1) - /ii/ (1) - /i 2 / (2) ) A 0i A (50)_idx 

+ / (g^ + g^ + g^ks - - k% 2 - a (,0)_i a { HO U dx 

+ J (p x (4fgW) - 7/i5 (3) - 2fgW) 6 1 A 2 A (h6)^dx 
+ J (D x (4fhW) - 7hhW - 2/^( 2 )) 0i A 2 A (<?0)_idx . 

This implies that every coefficient should be equal to zero, that is 

{ 2D x {f^f) - 5/i/(2) - 2//« + 2/^(3) + 2/^( 3 ) = 0; 

- J D 3 (2/< 7 ( 3 )) + Dl{2fgV> + 7/15^)) + D x (-2fgW _ 5 /i 5 ( 2 ) - 9/ 25 ( 3 ) - 2ghgW) 

+4/ 25 (2) + 2 ghgW + 5/ 35 {3) + 5(gh) l9 ^ + 3f l9 W - 9l fW - g 2 fW = 0; 
-L»3(2//t(3)) + D 2 x (2fhW + 7/^(3)) + D x (-2fhW - 5/^(2) _ 9/^(3) _ 2^(3)) 
+4/ 2 /i(2) + 2 5 /i/i(2) + 5/ 3 /i {3) + 5(#/i)i/i( 3 ) + 3/i/i( 1 ) - fci/W - /i 2 / { 2) = 0; 

+ 3(2)/^ + 5 (3) /i3 - /> {1) 5i - ^(2)52 - fc< 3 >ss = 0; 
4/5 {3) )-7/ l5 ( 3 )-2/ 5 (2)=0; 
, 4//i(3)) - 7/i/i(3) - 2//i(2) = . 

Substituting (40) into the above formulas and combining with (41), we obtain over-determined 
partial differential equations for the functions f(u\,u 2 ), Ph(u\,u 2 ), qh(ui, u 2 ), gh(uhU 2 ) and 
qh(ui,u 2 ). With the help of the package diffalg in Maple, we obtain the five cases listed in 
Theorem 4. o 



Acknowledgements 

The author would like to thank D.D. Holm, A.N.W. Hone, A.V. Mikhailov and P.J. Olver for 
inspiring discussions and useful comments. The author also thank the organizers of NEEDS 
2009 for the invitation and financial support, where the author began this research. 



References 

[1] Calogero F 1984 A solvable nonlinear wave equation. Stud. Appl. Math., 70:189-199. 



19 



[2] Dorfman I 1993 Dirac structures and integrability of nonlinear evolution equations. (Chich- 
ester: John Wiley & Sons Ltd.). 

[3] Fuchssteiner B 1983 Mastersymmetries, higher order time-dependent symmetries and con- 
served densities of nonlinear evolution equations. Progr. Theoret. Phys., 70(6):1508-1522. 

[4] Gel'fand I.M. and Dorfman I.Y. 1979 Hamiltonian operators and algebraic structures related 
to them. Functional Analysis and its Applications, 13(4):248-262. 

[5] Golenia J., Pavlov M.V., Popowicz Z., and Prykarpatsky A.K. 2010 On a nonlocal 
Ostrovsky-Whitham type dynamical system, its Riemann type inhomogeneous regulariza- 
tions and their integrability. Symmetry, Integrability and Geometry: Methods and Applica- 
tions, 6:002. 

[6] Holm D.D. 2006 Peakons. In J. P. Francoise, G.L. Naber, and S.T. Tsouv, editors, Ency- 
clopedia of Mathematical Physics, volume 4, 12-20 (Oxford: Elsevier). 

[7] Hunter J.K. and Saxton R.A. 1991 Dynamics of director fields. SIAM J. Appl. Math., 
51:1498-1521. 

[8] Hunter J.K. and Zheng Y 1994 On a completely integrable nonlinear hyperbolic variational 
equation. Physica D, 79:361-386. 

[9] Khesin B. and Misiolek G. 2003 Euler equations on homogeneous spaces and Virasoro 
orbits. Advances in Mathematics, 176:116-144. 

[10] Maltsev A.Y. and Novikov S.P. 2001 On the local systems Hamiltonian in the weakly 
nonlocal Poisson brackets. Physica D: Nonlinear Phenomena, 156(l-2):53-80. 

[11] Olver P.J. 1977 Evolution equations possessing infinitely many symmetries. Journal of 
Mathematical Physics, 18(6):1212-1215. 

[12] Olver P.J. 1993 Applications of Lie groups to differential equations, volume 107 of Graduate 
Texts in Mathematics, second edition. (New York: Springer- Verlag) . 

[13] Sakovich S. 2009 On a Whitham-type equation. Symmetry, Integrability and Geometry: 
Methods and Applications, 5:101. 

[14] Sergyeyev A. 2005 Why nonlocal recursion operators produce local symmetries: new results 
and applications. Journal of Physics A: Mathematical and General, 38:3397-3407. 

[15] Sanders J. A. and Wang J. P. 2001 Integrable systems and their recursion operators. Non- 
linear Anal, 47:5213-5240. 

[16] Sanders J.A. and Wang J.P 2001 On recursion operators. Phys. D, 149(1-2):1-10. 

[17] Wang J.P. 1998 Symmetries and Conservation Laws of Evolution Equations. PhD thesis, 
Vrije Universiteit /Thomas Stieltjes Institute, Amsterdam. 

[18] Wang J.P. 2009 Lenard scheme for two-dimensional periodic volterra chain. J. Math. Phys., 
50:023506. 



20 



